Electronic structure of periodic minimal surfaces 
— 'topological band structure' 
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Electronic band structure for electrons bound on periodic minimal surfaces is differential- 
geometrically formulated and numerically calculated. We focus on minimal surfaces because they 
are not only mathematically elegant (with the surface characterized completely in terms of "navels" ) 
but represent the topology of real systems such as zeolites and negative-curvature fullerene. The 
band structure turns out to be primarily determined by the topology of the surface, i.e., how the 
wavefunction interferes on a multiply-connected surface, so that the bands are little affected by the 
way in which we confine the electrons on the surface (thin-slab limit or zero thickness from the out- 
set). Another curiosity is that different minimal surfaces connected by the Bonnet transformation 
(such as Schwarz's P- and D-surfaces) possess one-to-one correspondence in their band energies at 
Brillouin zone boundaries. 
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I. INTRODUCTION 

There is a long history for the fascination in parti- 
cles bound on curved surfaces, which dates back to the 
early days of quantum mechanical There are two com- 
plimentary points of interests: one is how the particle 
motions are affected by the local curvature of the sur- 
face. The other is how the global topology (i.e., how 
the surface is wound) affects quantum mechanical wave 
functions. The latter problem becomes especially inter- 
esting if we consider a periodic surface embedded in the 
three-dimensional space. 

Geometrically, Schwarz, back in the 19th century, 
showed that we can make curved surfaces extend over the 
entire three-dimensional space by connecting hyperbolic 
(i.e., everywhere negatively curved) patches. Specifi- 
cally, Schwarz has constructed periodic minimal surfaces, 
where minimal means that the negatively-curved surface 
has a minimized area with the mean curvature (^(ki+k^) 
with Ki , K2 being the principal curvatures) vanishing ev- 
erywhere on the surface. B'B 

There are reasons from both mathematics and 
condensed-matter physics why periodic surfaces are in- 
triguing. First of all, periodic surfaces are of general 
interest in the condensed-matter physics, (i) In general, 
"crystals" (periodic structures) composed of surfaces are 
conceptually interesting as a class jpf periodic system on 
which electrons move. MackayQlJ has classified them 
group-theoretically (just as the ordinary crystals com- 
posed of atoms are classified with the space group), which 
he named 'flexi-crystallography'. 

(ii) In terms of materials science, periodic minimal sur- 
faces represent the topology of real condensed-matter sys- 
tems. These include not only conventional materials such 
as zeolites or a silica polymorph called melanophlogitetrEl, 
or isostructural silicon clathratesEEl, but recent advances 
in fabrication of exotic materials such as fullerenes or 



nanotubes have inspired further possibilities -such as 
negative-curvature fullerene (or Ceo zeolite)BEirb£l whose 
fabrication has been attempted with a zeolite as a 
templatelij. Their structures can be modeled as curved 
surfaces if we smear out atoms into a surface in the 
effective-mass sense, and it is a fundamental question to 
consider how a mobile (e.g., n) electron behaves on such 
surfaces 

Second, there are mathematical interests and simpli- 
fications when a periodic surface is minimal: we can 
exploit the Weierstrass representation, which enables us 
to specify the surface in a surprisingly simple manner 
in terms of "navels" . The representation also simplifies 
Schrodinger's equation as we shall show in the present 
paper. 

There are further mathematical fascinations specific 
to surfaces. One virtue of the structure constructed from 
surfaces is we can deform it. One can in fact deform one 
minimal surface into another with a differential geomet- 
rical transformation called the Bonnet transformation. 
We can then raise a question of how the band structure 
for one surface could be related with that for the trans- 
formed one. Another interest is that some-periodic min- 
imal surfaces, such as Schwarz's P-surfacal3, have a high 
symmetry ('interior-exterior' symmetry) that divides the 
space into two equivalent parts, which should be reflected 
in the electronic band structure. 

So in the present paper we address ourselves a ques- 
tion: how the electronic band structures should look like 
for periodic minimal surfaces. To start with, however, we 
have to envisage there are in general two ways (Figjlj) to 
make electrons confined to a surface: (a) One is to con- 
sider electrons bound to a thin_-curved slab of thickness d, 
where the limit d — > is takerO. (b) The other is to con- 
sider the surface with the degree of freedom normal to the 
surface ignored from the outset, i.e., a^-two-dimensional 
sheet is rolled into the curved surface&a. Either way it 
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has been shown that an effective potential arises from the 
curvature of the surface, but that the potential is differ- 
ent_between the two cases. Namely, the thin slab case 
(a)EZl has a potential 

-(fi 2 /8m)( Kl - K2 ) 2 , 

while the case (b)S has 

+ (H 2 /8m)(K 1 + k 2 ) 2 . 

The origin of the discrepancy was subsequently revealed 
by Nagaoka and coworkers xB when the degree of freedom 
normal to the surface is ignored, Dirac's prescription for 
constrained systems can be applied, but there is a room 
for ambiguity in the order of operators. If we adopt the 
conservation constraint, the resultant equation reduces 
to that in the d — > approach. When the surface is min- 
imal («i + «2 =0), the curvature potential is nonzero in 
general (since k% — k 2 = 2k % 7^ 0) in case (a), while the 
curvature potential vanishes identically in (b). 

For condensed-matter systems such as atoms arrayed 
along a curved surface, we should take the d — > ap- 
proach. Still, the difference in the band structure be- 
tween the two cases is curious. Namely, although we have 
a periodic system in either case, the periodicity imposed 
in case (a) is the periodicity in the strong potential that 
confines the electron into a thin slab (Fig. 1(a)), while in 
the case (b) the electron moves freely along the surface, 
where the only constraint is that an electron has to move 
in a space having a nontrivial topology. The topology 
can have a profound effect on the electron's wave func- 
tion, since, if we regard the periodic surface as a network 
of pipes (a cubic network for P-surface, diamond for D- 
surface, etc), the wavefunction interferes with itself along 
various paths wound around the "necks" . Thus the peri- 
odicity felt by an electron amounts to the strong confin- 
ing potential in case (a), while the periodicity only enters 
as a way in which the wave function interferes in case (b), 
and it is an intriguing question whether or not their band 
structures are similar. 

The purpose of this paper is, (i) to explicitly write 
down Schrodinger's equation for electrons on periodic 
minimal surfaces by exploiting Weierstrass's representa- 
tion in order to obtain the electronic band structure; (ii) 
to calculate and compare the band structures in cases (a) 
and (b). Unexpectedly the band structures turns out to 
be similar between cases (a) and (b), i.e., the bands are 
primarily determined from the topological way in which 
the wave function interferes with itself. The energy scale 
of the band structure (band splitting, such as a split of 
the 'd' band into E s and T2 ff , and band widths) is also 
universally ~ h 2 /2mL 2 with L being the linear dimen- 
sion of the unit cell of the periodic surface, (iii) We go 
further to 'Martensitic-deform' a surface to another con- 
nected by the Bonnet transformation. We shall show 
that there exists a curious one-to-one correspondence in 
their band structures, which illustrates another curious 
feature in the topological band structure. 



II. WEIERSTRASS REPRESENTATION OF 
MINIMAL SURFACES 



We start with a mathematical prerequisite for rep- 
resenting minimal surfaces. A two-dimensional surface 
r(g 1 ,g 2 ) embedded in a three-dimensional space can be 
expressed in terms of two dimensional coordinates g 1 , g 2 , 
where (g x ,g 2 ) = (u, v) are called isothermal when the 
metric tensor g^ is diagonal with 

dr ■ dr = gn(dudu + dvdv). 

What Weierstrass and Enneper have found is that a nec- 
essary and sufficient condition for r(u, v) representing a 
minimal surface with isothermal (u, v) € S (S: a simply 
connected region) is that there exist F, G, functions of 
w = u + iv, with which r(u, v) = (x(u, v), y(u, v), z(u, v)) 
is expressed as 



r(u,v)=Re(f F(l~G 2 )dw 
/ iF(l + G 2 )dw, / 2FG dw 



(1) 



where Wq is a constant, and FG 2 is assumed to be regular 
(i.e., m-th-goles of G assumed to coincide with 2m-th ze- 
ros of F) .□'B'Ej If there are singularities that violate this 
condition, we can exclude these points by incising cut(s) 
to make S a Riemann surface. Thus there is a one-to- 
one correspondence between a minimal surface and the 
functional form of F, G. 

Now, Schrodinger's equation for a curved surface, ex- 
pressed with two dimensional coordinates (g\g 2 ) and 
metric tensor g.-,-, is written as 
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where summations over repeated indices are assumed. 
This equation is for model (a), while we can replace 
the second term in the bracket (potential term) by 
( Kl + K2 ) 2 for model (b). 
In the Weierstrass-Enneper representation, every 
quantity in the Schrodinger's equation can be expressed 
in terms of F and G, since the Laplacian, the first term 
in the angular brackets in eq.|^, reduces to (d 2 /dq 2 )/ ' yfg 
in the isothermal coordinates, where 



g = det{.g.y} = 



1 



|F|(|G| 2 



while we can plug in k x = -k 2 = 4|G'|/|F|(|G| 2 + l) 2 for 
the curvature term. Schrodinger's equation for periodic 
minimal surfaces then reduces to 
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where e = E/ (h 2 /2m). As evident from the Weierstrass 
representation (|l]), F has the dimension of length and G 
is dimensionless. Hence the energies in minimal surfaces 
always scale as E/(h 2 /2mL 2 ), where L ~ F ~ linear 
dimension of the unit cell (a precise expression given in 
cqn.|| below). 

We can have a more transparent form when G(w) = w 
(as is often the case with periodic minimal surfaces, in- 
cluding P-surface). In this case we can exploit the stereo- 
graphic map (Gauss map) from the infinite complex plane 
(u, v) to a unit sphere (6, (j>), 
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After a bit of algebra, we finally arrive at the differential 
equation for (8, (ft), 
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Curiously, a common coefficient (1 — cos6>) 4 /\F\ 2 fac- 
tors out for the Laplacian (the first three terms in the 
large parentheses in the above equation) and the curva- 
ture potential (+1), which is made manifest due to the 
Weierstrass representation. Hence the kinetic and poten- 
tial energies vary in a correlated manner as we go from 
one minimal surface to another by changing F . However, 
this does not imply that the curvature potential always 
exerts as large an effect as the kinetic energy does, since 
the expectation values of the kinetic and potential ener- 
gies depend on the amplitude of the wave function, and 
a quantitative study is required. 



III. RESULTS FOR SCHWARZ'S P-SURFACE 

So we first take Schwarz's P-surface as a typical triply- 
periodic (periodic in x, y, z) minimal surface, or a 'sim- 
ple cubic' minimal surface (belonging to a space group 
Im3m inJ.he.language of the flexi-crystallography) . Some 
authorscHt 2 ] gave the Weierstrass-Enneper representa- 
tion for the P-surface as 



F{w) = iL/y/1 - 14w 4 



(5) 



where L ~ linear dimension of the unit cell (to be precise 
the unit cell size is 2.157L, which is given as an elliptic 
integral) . 

A unit cell of the P-surface comprises eight identical 
patches, as shown in Fig.||. With the stereographic map- 
ping discussed above, a unit cell is mapped onto two 
spheres, connected into a Riemann surface via four cuts, 



which have to be introduced to make FG 2 non-singular 
since F has poles. 

Differential geometrically F is in general specified 
solely by such poles in a form F oc 1/ Jli( w— w i) n where r\ 
determines the topology of the surface. In fact the poles 
correspond, in the language of differential geometry, to 
the navels (umbilical points), which are defined as the 
points where every cross section is inflected, with the two 
principal curvatures becoming degenerate, k\ = K2 (= 
for a minimal surface). So a periodic minimal surface is 
completely characterized by the navels that appear peri- 
odically. The curvature potential cx — (aci — K2) 2 hi case 
(a) varies on the surface, which may be called a 'crystal 
field' in flexi-crystals. Navels then specify the positions 
where the curvature potential becomes maxima (= 0), 
while the minima occurs at the maxima of the absolute 
value of Ki(= —K2 for a minimal surface). In the P- 
surface the navels (potential maxima) occur at eight 'Af- 
fensattel' (monkey's-saddle) points in a unit celled, while 
the potential minima occur at four poiute around each 
nape of the neck as depicted in Fig.^(a) .cJ 



IV. RESULTS FOR THE BAND STRUCTURE 

Band structure for case (a) In Schrodinger's equa- 
tion for periodic minimal surfaces, eqn.(|j), the variables 
6, 4> cannot be separated, so that we have solved the equa- 
tion numerically by discretizing 6, <fi to diagonalize the 
Hamiltonian matrix. In discretizing the spherical coordi- 
nates, a special care is taken around the navels, since the 
Jacobian J of the transformation to the Gauss sphere is 
singular there (Fig.||) . The band structure is obtained by 
connecting the adjacent unit cells with appropriate phase 
factors. 

We now come to the result for the band structure for 
the P-surface in Fig.^j (curves), and typical wavefunctions 
at r and H in the bcc Brillouin zone (r in the simple- 
cubic zone) in Fig.|^. The P-surface happens to divide 
the space into two equivalent parts, since a body center 
enclosed by the surrounding unit cells has the same shape 
as the original unit, so that we first note that the true 
symmetry is body-centered cubic rather than simple cu- 
bic. So the bands are displayed on the Brillouin zone for 
bcc. In accord with the above argument, the energy scale 
(band width, splitting etc) is ~ Ti 2 /2mL 2 . This is of the 
order of 1 eV for L ~ 10 A, the unit cell sizje-assumed for 
a hypothetical negative-curvature fullerenco. 

The curvature, or the effective mass, of these bands 
are either positive (electron-like) or negative (hole-like) 
according to the nature of the wave function. The mass 
cannot be estimated with a simple k • p perturbation, 
since the perturbation cx k • p derives from the fact 
that Hq oc p 2 while Hq has no such simple form on a 
curved surface. In other words, the k • p formula has 

Y l i{{ i \pA3){j\Pv\ i ))l{ E i - E i)i so we would have to cal- 
culate the matrix elements of p for wavefunctions that 
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are finite only along the surface. 

Band structure for case (b ) Now we are in position 
to compare the (a) confinement case (curves in Figj4|) 
with the (b) rolled case (dotted lines in Fig|5|). We can 
immediately see that the two band structures are rather 
similar up to some offset (2.34?i 2 /2mL 2 ). This is surpris- 
ing, since there is no apriori reason why they should be. 
To be more precise quantitative features characterizing 
the band, i.e., the effective mass and band widths, are 
similar between the two cases. So we conclude that the 
band structure does not essentially depend on the way in 
which electrons are confined, at least for 'gently' curved 
surfaces such as minimal surfaces where the there are no 
sharp edges that would give large curvature potentials. 



V. BONNET TRANSFORMATION 

The next important question is: are the band struc- 
tures for surfaces connected by the Bonnet transforma- 
tion related? 

Bonnet transformation The deformation of the P- 
surfacc to other periodic minimal surfaces can be imple- 
mented by the Bonnet transformation, which is confor- 
mal and is represented by an elliptic transformation. A 
beautiful asset of the Weierstrass representation for mini- 
mal surfaces is that the Bonnet transformation is simply 
represented by a phase factor, F — > Fe 1 ^ , in eqn.(|), 
where is called the Bonnet angle. If we apply ptius 
to the P-surface (cubic), the transformation changesEJ it 
into the G-surface (gyroid with (3 = 0.21 In) and the D- 
surface (diamond with (3 = n/2), which may be regarded 
as a 'Martensitic transformation' in the words of Rcf.EZl. 
The structure of the D-surface is depicted, along with the 
P-surface, in Fig.|| Since the D-surface has a diamond 
symmetry, its unit cell contains two 'cages'. 

We can first note that the Bonnet transformation 
preserves the metric tensor and the Gaussian curva- 
ture. This implies that all the surfaces connected by the 
Bonnet transformation obey the identical Schrodinger 
equation within a unit patch. Indeed, if we look at 
Schrodinger's eqn.(^), F only enters as so that 

F — > Fe 1 ^ does not alter the equation. Although this is 
curious enough, this does not mean that the band struc- 
tures are identical, since, while the transformed surfaces 
share the same genus (= three for P and D), the way in 
which unit cells are connected is different among them. 

Band structure of a Bonnet-transformed surface 

Figure || compares the band structures for the P- and 
D-surfaces. The two band structures are indeed differ- 
ent due to the difference in three dimensional connection 
of the unit cells discussed above. Curiously, however, 
we find that the values of the band energy at special 
points (Brillouin zone corners, edges and face-centers) 
have identical set of values between different surfaces. 
Namely, a close comparison of the two band structures 
reveals that the band energies exactly coincide, where the 



'law of correspondence' is 
P-surface D-surface 

r, h & tTr 

N ^ X, M 
This can be explained from the property of the Bonnet 
transformation that does not change Schrodinger's equa- 
tion. For this purpose we have to look at the unit cells 
more closely. In Fig.^(a), we show how unit patches are 
connected for the D- and P-surfaces. We have indicated 
in the figure how the eight patches in a unit cell, num- 
bered with 1 through 8, are connected to other patches in 
the adjacent cells by marking the edges with those num- 
bers (i.e., if an edge is marked with, say, 7, the adjacent 
patch should be 7). The wave function should be con- 
tinued to the adjacent patch with a certain phase factor. 
We have indicated the connection coefficients, 

Pi = exp(j^), 

where <pi is the Bloch phase along i = x,y, z. Hence this 
diagram fully characterizes how the Bloch wave functions 
are connected on the periodic minimal surfaces in terms 
of patches. We can then compare the coefficients for the 
P- and D-surfaces to extract a correspondence at special 
points in fc-space. 

The diagrams introduced here have naturally differ- 
ent patch numbers assigned on them between the P- and 
D-surfaces, since the way in which the patches are con- 
nected (i.e., numbers attached to the edges) is different 
between them. However, we can make them identical, 
if we rearrange the connection numbers by noting the 
symmetry. Since the parity inversion o~ p with respect to 
the center of a unit cell preserves the Hamiltonian and 
does not change the special fc-vectors on the zone edges 
either, the group theory dictates that the eigenstates on 
those fc-points should have the parity 1 or -1. Then we 
can always construct the wave function for one half of 
the unit cell by multiplying the wave function for the 
other with 1 or -1. Namely, if we employ a simple cu- 
bic unit cell for the P-surface (Fig.*, which is twice the 
bcc unit cell depicted in Fig.||(a)) to make the correspon- 
dence clearer, an application of a p to the upper half?, the 
lower turns out to have just the same connection num- 
ber s? as in the D-surface as shown? in Fig.||(b). So 
we have now the same connection numbers between P 
and D, where the only difference is different connection 
coefficients (pi's) between P and D as indicated in the 
figure. If we compare these, we end up with a 'law of 
corresponding fc-points', 

P-surface D-surface 
fc-points^, jo y ,jo z ) k-pomts(p x ,p y ,p z ) 



r(W) 



i 

-l & 



F(l,l,l) 
R(-l,-l,-l) 



M(-l,-l,l) op = -1 o X(l,l,-1) 
<T p = 1 O M(-l,-l,l) 
We can immediately translate this into the correspon- 
dence found above, r, H <=> T, R, etc, if we note the 
relation, 
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simple cubic bcc 
T <5 r, H 
M ^ N 
between the simple-cubic and bcc unit cells. 

The wave functions shown in Figs.0, are actually re- 
lated through this relation. A simplest way to confirm 
this is to note that the wave function on each unit patch 
behaves in a similar manner. In fact, the wave function 
4>) hi eqn.[| is identical between the two surfaces. 

VI. DISCUSSIONS 

The band structures revealed here should have impor- 
tant implications on various physical properties. These 
should include transport properties as well as the cy- 
clotron resonance, which can detect the effective mass 
arising from topological band structures. Since the mass 
is determined by the interference of wave functions, ef- 
fects of external magnetic fields should also be interest- 
ing. 

We can finally comment that, i£ we adopt foams of 
graphite to realize curved surfacesQ, then the equation 
of motion of 7r electrons on the network of the honey- 
comb lattice will become, in the effective mass picture, 
the problem of zero-mass Dirac equation (i.e., Weyl's 
equation) on curved surfaces. While we have ignored 
spin degrees of freedom here, the spin connection on the 
surface will give rise to a Berry's geometrical phase. 

H.A. wishes to thank Alan Mackay for illuminating cor- 
respondences, Hiroshi Kuratsuji and Koichi Fukuda feu- 
discussions, and Yasuo Nozue for pointing out the ref.E-3. 

FIG. 1. Two ways to prepare a surface: (a) to introduce po- 
tential barriers that confine electrons in a three-dimensional 
space to a thin membrane, or (b) to roll a sheet of free elec- 
trons into the curved surface. 

FIG. 2. The structures of P(a) and D(b) surfaces. We 
show a unit patch on the left panel, and a full unit cell on the 
right. The grey-scale in the right panel represents the curva- 
ture potential, where shaded (open) circles depict potential 
minima (maxima, coinciding with the navels), (c) The stere- 
ographic projection from minimal surface to Gauss spheres. 
A unit patch of P or D corresponds to a pair of 1/8 spheres. 

FIG. 3. Descretization for the spherical coordinates: (a) 
a mesh with even intervals in 6 and <j> and (b) uneven ones, 
adopted here, which take care of the singular point in the 
Jacobian. 

FIG. 4. The energy band structure in units of Ti 2 /2mL 2 
is shown for the P-surface, when the curvature potential is 
considered (curves) or ignored (dots) with an energy offset to 
make the band bottoms coincide between the two cases. The 
inset depicts the Brillouin zone for a bcc unit cell. 



FIG. 5. Typical wavefunctions for a unit cell of P-surfaces 
with positive (negative) amplitudes color-coded in red (blue) . 
Their eigenenergies are indicated in Fig.^. 

FIG. 6. The band structures for the P- and D-surfaces. 
Horizontal lines indicate how the energies at the zone center 
or edges coincide between the two cases. 

FIG. 7. Typical wavefunctions for a unit cell of the Bon- 
net-connected P- and D-surfaces with positive (negative) am- 
plitudes color-coded in red (blue). Their eigenenergies are 
indicated in FigQ 

FIG. 8. (a)The way in which patches (labeled by large 
numbers) are connected to those in adjacent unit cells are 
indicated by small numbers attached to the edges for the 
D-surface. The Bloch phase factors (pi's) are also shown. In 
the left panels the patches are flattened and expanded, while 
the right panels depict the actual three-dimensional shapes. 
(b)The corresponding diagram for the P-surface, where we 
have rearranged the numbers to make them identical with 
those for the D-surface by exploiting the symmetry. Accord- 
ingly the Bloch phase factors for the P-surface involve the 
parity inversion (<r p ). An example is shown in (c), which in- 
dicates how patches 7,8 are neighboring 3,4 through p x a p . 
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